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Abstract: In the present work, we use the homotopy perturbation method (HPM) to solve the Newell- 
Whitehead-Segel non-linear differential equations. Four case study problems of Newell- Whitehead- 
Segel are solved by the HPM and the exact solutions are obtained. The trend of the rapid convergence 
of the sequences constmcted by the method toward the exact solution is shown numerically. As a result 
the rapid convergence towards the exact solutions of HPM indicates that, using the HPM to solve the 
Newell-Whitehead-Segel non-linear differential equations, a reasonable less amount of computational 
work with acceptable accuracy may be sufficient. Moreover the application of the HPM proves that the 
method is an effective and simple tool for solving the Newell-Whitehead-Segel non-linear differential 
equations. 
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INTRODUCTION 

Recently lots of attentions are devoted toward the semi-analytical solution of real-life mathematical 
modeling that is inherently nonlinear differential equations. Most of the nonlinear differential equations do not 
have an analytical solution. The idea of the homotopy perturbation method was first pioneered by (He, 1999) 
Later the homotopy perturbation method (HPM) which is a semi-analytical method is applied to solve the non- 
linear non-homogeneous partial differential equations (He, 2005; Yildirim, 2009; He, 2006; He, 2006; He, 2006; 
He, 2005; He, 1998; He, 1998; He, 1999; He, 1999; He, 2000; He, 2006; He, 2005; He, 2005; He, 2005; Ko<?ak 
et al., 2011; Gepreel, 2011; Cao and Bo. Han, 2011; Wazwaz, 2009). Ezzati and Shakibi (Ezzati, 2011) solved 
the Newell- Whitehead equation using the Adomian decomposition and multi-quadric quasi-interpolation 
methods. They concluded that the Adomian decomposition and multi-quadric quasi-interpolation methods are 
reasonable methods to solve the Newell- Whitehead equation with acceptable accuracy. 

In the present work, the homotopy perturbation method (HPM) is applied to obtain the closed form 
solution of the non-linear Newell-Whitehead-Segel equation. Four case study problems of non-linear Newell- 
Whitehead-Segel equations are solved using the HPM. The trend of the rapid convergence towards the exact 
solution is shown when compared to the exact solution. The Newell-Whitehead-Segel equation models the 
interaction of the effect of the diffusion term with the nonlinear effect of the reaction term. The Newell- 
Whitehead-Segel equation is written as: 
c? 2 ii 

— = k — - + au - bu q (1) 

Ot OX* 

Where a? b and k are real numbers with k > and q is a positive integer. In Eq. (1) the first term on the 

&u y 

left hand side, — , expresses the variations of with time at a fixed location, the first term on the right 

hand side, — - expresses the variations of with spatial variable x at a specific time and the remaining 

terms on the right hand side, ati — , takes into account the effect of the source term. In Eq. (1) t) is 

a function of the spatial variable x and the temporal variable t, with X ■ ER and t >0. The function u(:r,t) 

may be thought of as the (nonlinear) distribution of temperature in an infinitely thin and long rod or as the flow 
velocity of a fluid in an infinitely long pipe with small diameter. The Newell-Whitehead-Segel equations have 
wide applicability in mechanical and chemical engineering, ecology, biology and bio-engineering. 
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2. The Idea Of Homotopy Perturbation Method: 

The homotopy perturbation method (HPM) is originally initiated by (He, 1999; He, 2005; Yildirim, 2009; 
He, 2006; He, 2006; He, 2006; He, 2005; He, 1998; He, 1998; He, 1999; He, 1999; He, 2000; He, 2006; He, 
2005; He, 2005; He, 2005; Kogak et al., 2011). This is a combination of the classical perturbation technique and 
homotopy techinique. The basic idea of the HPM for solving nonlinear differential equations is as follow; 
consider the following differential equation: 

= Oj- (2) 

Where E (it) is any differential operator. We 

H(u,p) = (1 -p)(F(tt) - if (it)) + p(ff(u.)) = 0. construct a homotopy as follow: 



(3) 

Where are F(u), H(u) are functional operators with the known solution v 0 . It is clear that when p is equal to 
zero then H(u, 0) = F(u) - H(u) = 0, and when p is equal to 1 , H(u, 1) = E(u) = 0. It is worth nothing that is the 
embedding parameter p increases monotonically from zero to unity the zero order solution v 0 continuously 
deforms into the original problem E(u) = 0. The embedding parameter, p £ [0,1], is considered as an expending 
parameter(He, 1999). In the homotopy perturbation method the embedding parameter p is used to get series 
expansion for solution as: 

DO 

u = '^p , v i = v 0 + pv 1 +p 2 v 2 +p 3 v 2 -i (4) 

i =0 



When p — M, then Eq. (3) becomes the approximate solution to Eq. (2) as: 
u = v a + v 1 + v 2 + v 3 H (5) 

The series Eq. (5) is a convergent series and the rate of convergence depends on the nature of Eq. (2) (He, 
1999; He, 2005; Yildirim, 2009; He, 2006; He, 2006; He, 2006; He, 2005; He, 1998; He, 1998; He, 1999; He, 
1999; He, 2000; He, 2006; He, 2005; He, 2005; He, 2005; Kogak et al., 2011). It is also assumed that Eq. (3) has 
a unique solution and by comparing the like powers of p the solution of various orders is obtained. These 
solutions are obtained using the Maple package. 

3 . The Newell-Whitehead-Segel Equation : 

To illustrate the capability and reliability of the method, four cases of nonlinear diffusion equations are 
presented. 



Case I: 

In Eq. (1) for a = 2 f b = 3 x k = 1 and q = 2 the Newell-Whitehead-Segel equation is written as: 



at 



d 1 



U 



— = — — - + 2u - 3 u 2 



dx 2 

Subject to a constant initial condition 

u(x,0) = X, 



( 6 ) 

(7) 



We construct a homotopy for Eq. (6) in the following form: 



— r-2u + 3v 2 



H(v,p) = (1 - p) — - 



~dv 


du 0 - 


+ p 


dv 


d 2 v 




dt . 


dt 


dx 2 



= 0 



( 3 ) 



The solution of Eq. (6) can be written as a power series in p as: 
v=v a +pv 1 + p 2 v 2 +■■■ 



(9) 
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Substituting Eq. (9) and Eq. (7) into Eq. (8) and equating the terms with identical powers ofp: 





dv a 


dUfy 




£ 

O 

X 

O 

II 




p a : 


dt 


dt ’ 




A, 




5 'ii. 


du n d 2 v n 








P 1 


± 

dt 


dt dx 2 


1^0 , 


£ 

H- 

X 

O 

II 


0, 




dv 7 


d 2 v 1 


+ u 1 (2 - 3u 0 ) , 


u 2 (x,0) - 




p 2 : 


£ 

dt 


= dx 2 3 "0"1 


0, 






d 2 v 2 


— 3u 2 +u 2 (2 — 3u 0 ) , 






p 2 : 


d 

1 ~dt 


= dx 1 3 “» U2 


^3 (t,0) = 


0. 




Using 


; the Maple package to solve recursive sequences, Eq. 


(10), we obtain the followings: 



= A, 

i.'j (x, t) = A(2 — 3 A) t, 

v 2 (x, t) = 21(2 — 31) (1 — 31) ^ , (11) 

t 3 

u 3 ( x , t) = 21(2 - 31)(271 2 - 1S1 + 2) — . 

By setting p = 1 in Eq. (9) the solution of Eq. (6) can be obtained as 
u = u 0 + v 1 +■ v 2 +■ +■ ... Therefore the solution of Eq. (6) is written as: 

v (x ,t) = 1 + 1(2 - 31) t + 21(2 - 31) (1- 31) £ + 21(2 - 3l)(27l z - 181 + 2) ^ + 



The Taylor series expansion for 






--+A-Ae at 



is written as: 



-As 



2t 



■ I - 1 — As 



it 



= 1 + 1(2 - 31) t + 21(2 - 31) (1 - 31) 



2! 



+ 21(2 - 31) (271 2 -181 + 2) — + ■■■ 

By substituting Eq. (13) into Eq. (12), thus Eq. (12) can be rewritten as: 
-^le 2t 

v (x r t ) = - 

--+ A- Ae 2t 



( 12 ) 



(13) 



(14) 



This is the exact solution of the problem, Eq.(6). Table 1 shows the trend of rapid convergence of the 
results of 5 0 (^ t) = (x, t) to 5 5 (x f t) = using the HPM. The rapid convergence of the 

solution toward the exact solution, the maximum relative error of less than 0.0046%, is achieved as shown in 
table 1. 

Table 1 shows: he percentage of relative errors of the results of S 0 (r # t) = v 0 (r # t) to 
5 5 (x 0 = Sp=oVi Of, fc) of the HPM solution of Eq. (6) for X=0.1 

Case II: In Eq. (1) for a — 1, b =■ 1, k = 1 and (7 — 2 the Newell-Whitehead-Segel equation is 
written as: 
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Table 1: 



percentage of relative error (%RE) 


t = 0.1 


S a {x,t) 


0.1540788596 




Siix, 0 


0.01027226573 




S 2 (x,t) 


0.0002058041542 




Si (x,t) 


0.00001949256084 




S 4 (x r t) 


0.000001708760704 




S 5 {x,t) 


3.214500334 e-8 


t = .3 


S 0 (.x,t) 


0.3835101093 




S^x, t) 


0.06910026503 




S 2 (x..t) 


0.003074197736 




S 3 (x,t) 


0.001358981068 




S 4 (_x, t) 


0.0003091665977 




S B (x,t) 


0.00001230390524 


II 


S a (x,t) 


0.4680703804 




S ± (x, 0 


0.1063582391 




S 2 {x,t) 


0.005078839539 




S 3 (x,t) 


0.003988077959 




s 4 (x, 0 


0.001125247109 




S s ,(x,t) 


0.00004585924829 



du 

dt 



d 2 u 

= JP + U ~ U 

Subject to initial condition 



w(x,0) = 




To solve Eq. (15) we construct a homotopy in the following form: 



H(v t p) = (1-p) 



"3d 


du 0 - 


+ p 


dv 


d 2 v 


.dt 


dt . 


dt 


dx 2 



~ — 7 — u(l - u) 



= 0. 



( 15 ) 

( 16 ) 



( 17 ) 



The solution of Eq. (15) can be written as a power series in p as: 
v =v a +pv 1 + p 2 v 2 +-- (IS) 

Substituting Eq. (18) and Eq. (16) in to Eq. (17) and equating the term with identical powers of p , leads to 

1 

u °(*'°) = t TT7 

(l + e^\ 



P Q -. 


d "2-Iq 5’W- q 

dt dt’ 


■H 

5 


dv 1 du Q 


P ■ 


dt dt 




dv 2 


p . 


dt dx 2 


p 3 = 


3li 3 d 2 v 2 


dt dx 2 



(19) 
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Using the Maple package to solve recursive sequences, Eq. (19), we obtain the followings: 









u z (x,t) = — 




(20) 



v 3 (x, t) = 



By setting p — 1 in Eq. (18) the solution of Eq. (15) can be obtained as 
v = + v ± + v 2 + ^3 + — Therefore the solution of Eq. (15) is written as: 



/ 3/ -fV B 

f 1 + 8 % * j f 1 +B'* j 




The Taylor series expansion for 






s is written as 



(i + .**) (i + .*) 3 (i + .*J 

125 ( 4 ( e ^) 2 - 



7e^ + 1 



216 



(l + ^) B 




( 21 ) 



( 22 ) 



By substituting Eq. (22) into Eq. (21), the Eq. (21) can be reduced to 
vQx, t) = 1 



(»M) 



(23) 



This is the exact solution of the problem, Eq. (15). Table 2 shows the trend of rapid convergence of the 
results of 5 0 (x r t) = Vq (:r , t) to = Zf= q ^ t} using the HPM solution toward the exact 
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solution. The maximum relative error of less than 0.0008% is achieved in comparison to the exact solution as 
shown in table 2. 



Table 2: 







percentage of relative error I 


(%RE) 






t=I 

II 

* 


x = l.S 


% = 1.8 


t = 0.1 


5 0 (x,t) 


0.09374691950 


0.1010110875 


0.1051585828 




Site, 0 


0.003020748923 


0.003848096613 


0.004359470571 




s 2 (x,t) 


0.00001102987604 


0.00002050070585 


0.00004357185889 




S 3 (x,t) 


0.000003008522247 


0.000003116144165 


0.000002949567767 




S 4 (x, t) 


3.247853526 e-8 


7.597179054 e-8 


1.007266564 e-7 




5 5 (x,t) 


2.475184340 e-9 


6.284562692 e-10 


4.103584300 e-10 


II 

U-j 


S 0 (x,t) 


0.2480807805 


0.2663112457 


0.2766522515 




5ife 0 


0.02225395729 


0.02841934500 


0.03220845585 




S 2 (x 1 t) 


0.0003852758255 


0.0003051199126 


0.0008095286291 




S 3 (x,t) 


0.0002055820118 


0.0002152872468 


0.0002058262157 




S 4 (x, t ) 


0.000005574634536 


0.00001431284116 


0.00001929380478 




S 5 0t,t) 


0.000001472642781 


0.000001037290051 


5.725167943 e-7 


r-t 

II 

4^ 




0.3115489693 


0.3338572016 


0.3464690907 




s ± (x, t) 


0.03586206901 


0.04586954916 


0.05200200785 




S 2 (x,t) 


0.0009882331330 


0.0004901161293 


0.001569416719 




S 3 (x,t) 


0.0005982451688 


0.0006298759390 


0.0006050514464 




S 4 (x, t) 


0.00001947940975 


0.00005317431566 


0.00007241822783 




S s {x,t) 


0.000007711072685 


0.000005555853236 


0.000003218048969 



Table 2 shows: the percentage of relative errors of the results of 5 0 ti, 0 = i? 0 U, fc) to 
5 5 (x t) = Gr, t) of the HPM solution of Eq. (15). 

Case III: In Eq. (1) for a = 1, b = l f k = 1 and q = 4 the Newell-Whitehead-Segel equation 



becomes: 
du d 2 u 
dt 



— = — ^ + u- u 



dx 2 

Subject to initial condition 



(24) 






BJf 

\1 + e-i™ 



2 




(25) 



We construct a homotopy for Eq. (24) in the following form: 



= (i-p) 



~dv 


du Q - 


+ P 


dv 


ft. 


[dt 


dt J 


dt 


1 

Cl, 

r-j 

1 



■V + v 4 



= 0. 



(26) 



The solution of Eq. (24) can be written as a power series in p as: 
v =v Q + pv 1 + p 2 v 2 + ■■■ 



(27) 



Substituting Eq. (27) and Eq. (25) into Eq. (26) and equating the terms with identical powers of p: 
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Q m 

^ dt dt 



U Q (*,0) = 



(l + e-^y 



1 dv ± d 2 Ufy du 



dt dx 

dv 2 e? 2 !^ 
dt dx 2 



3 _ 3u 3 d z v 2 



2 -^r+v bd- v o 3 X 



p 2 ^ 



dt 

^(J 3 ^ + V i _ , 






v 1 (v, 0} = 0., 

v 2 ( x , 0 ) = 0 , 
u 3 (x, 0 ) = 0 ., 



(28) 



Using the Maple package to solve recursive sequences, Eq. (28), we obtain the followings: 
a 

1 V 






'x'l + 






vA x >0 = g 



3 

J — — i 



\ 






(l + e--Uy 



/ 






u z (x, t) = 



49 

50 



s \ 

-=i 1 






(t + e^y j 



(29) 



-27e^*+9^ 



v 3 ( x r t) = 



343 

1000 






27e' JLCf + 9 



( i+e U‘ / 



By setting p — 1 in Eq. (27) the solution of Eq. (24) can be obtained as 
u = u 0 + v 1 + u 2 + + ... thus the solution of Eq. (24) can be written as: 






vCr,t)=( ± 3 -)+^ 

\1 + eJH 1 ) 5 



\ 



e -.'n> 



49 
t + — 

50 



+ 



343 

1000 



^26^* - 3 j 



^i + y y y ^i + e ^ y j 

^4^6’^®*^ — 27e'^ x + 9^e , > fi °*^ ^_ 3 

V (i+^) T j 



(30) 
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2 

The Taylor series expansion for ^ tank ^ ^ J + 



is written as: 




Comparing Eq. (31) with Eq. (30), thus Eq. (30) can be reduced to 



2 , 




(32) 



This is the exact solution of the problem, Eq. (24). Table 3 shows the trend of rapid convergence of the 
results of 5 0 (^t) = to 5 s (^t) = using the HPM solution toward the exact 

solution. The maximum relative error of less than 0.007% is achieved in comparison to the exact solution as 
shown in table 3. We can conclude that the HPM is one the most suitable and friendly method in solving the 
Newell-Whitehead-Segel equation. 

Table 3 shows: the percentage of relative errors of the results of 5.j(x t) to 

t) = Ef =0 i^ Gx t) of the HPM solution of Eq. (24). 



Case IV: 

In this case we will examine the Newell-Whitehead-Segel equation for a = 3> b = 4, k = 1, q = 3, 
chi d 2 u 
dt 



— = — — ^ + 3 u - 4 u , 



dx 2 

Subject to initial condition 



(33) 



We construct a homotopy for Eq. (33) in the following form: 

H(v r p) = (1-p) 

The solution of Eq. (33) can be written as a power series in p as: 
Li = V ± + p 2 V 2 H 



~dv 


du a - 


+ p 


dv 


d 2 v 


3t 


dt J 


dt " 


dx 2 



= 0 . 



(34) 

(35) 

(36) 



Substituting Eq. (36) and Eq. (34) into Eq. (35) and equating the terms with identical powers of p : 
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,dv 1= du 1 

p ■ st dt' 






A _ 



P 



d Vl d 2 v 0 dv Q 

■it 



(r, 0) = 0,, 



2 c?u 2 3 z u 1 


1 2 u 4- Tii 


u 2 < 

12d 2 D(, d 3 ( 


:x,0) 

x,0) 


= o„ 
= 0. 


^ dt dx 2 

3 . 9u 3 _ ^2 


LjL'Uq eJ-^ | J 1 L'^ 

-1- ^11 1 ?7I 


^ dt dx 2 

Table 3: 


T J L 1 -? JL^cJqcJ^ 


percentage of relative error (%RE) 




X = 1 


x = l.S 




x = 1.8 


t = 0.1 5 0 Ut) 


0.09323514238 


0.1045372618 




0.1099789256 


Sxix, t) 


0.001725301473 


0.003516790607 




0.004498961488 


S 2 {x,t) 


0.0002100221036 


0.0001216772460 




0.00005170881962 


S 3 (x,t) 


0.00001173612982 


0.00001539927439 




0.00001472293353 


S 4 (x, t) 


3.620778386 e-7 


2.929797495 e-7 




5.746720664 e-7 


s 5 0,0 


5.414213656 e-8 


3.614791446 e-8 




1.434500497 e-8 


t = .3 5 0 Otr, t) 


0.2396054415 


0.2702774140 




0.2852064054 


S t (_x, t) 


0.009390569965 


0.02330928345 




0.03106724554 


S 2 (x > t) 


0.005215736776 


0.003375969333 




0.001825357483 


S 3 (x,t) 


0.0007262150389 


0.001037577469 




0.001023346798 


S 4 ix, t) 


0.00009555513293 


0.00004044412369 




0.0001029640746 


S 5 (x,t) 


0.00003280557979 


0.00002473085989 




0.00001198739229 


II 

a- 


0.2962707966 


0.3351845099 




0.3542406943 


Si(x, t) 


0.01219208796 


0.03518330953 




0.04811460219 


S 2 (x,t) 


0.01183961035 


0.008037424274 




0.004713588643 


S 3 (x,t) 


0.001990820701 


0.002987593293 




0.002996129616 


S 4 (x r t) 


0.0004128329726 


0.0001164720814 




0.0003682018367 


5 5 (x,0 


0.0001681121669 


0.0001337455766 




0.00006941923143 



Using the Maple package to solve recursive sequences, Eq. (37), we obtain: 



u 0 O,t)= - — 

-Jsx 



-J&x 






VifoO =2 



(e V6je + e = *J 



t. 



, , 81 3 

u z (x,0 = — T 



e ^ e V 



+e»*) t z 



( 1±\ 

/ e V6* + g j 
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v 3 (x, t) = 



243 3 



a ' 






16 J 4 



|e V&ae + e a 



By setting p — 1 in Eq. (36) the solution of Eq. (33) can be obtained as 
v = v Q + v ± + v 2 + v 3 + ... 

Thus the solution of Eq. (33) can be written as, 

v ft, t) = 



i? . > I? , , 3i 

^ .T'+(.^)’ t(.f J 



243 

16 



(39) 



The Taylor series expansion for 



1 r- (fx-lt) 



is written as 



V fcJt’ 



_|_ g V 2, * 3 7 



Vfoe 



3 e 

4 _ _ 

gV&s -[- g *a * 



9 

+ T 



_ V t> 






+ 



31 +.¥“),» 
4 j4 . _ 

3 V&Je 



243 

TfT. 



^e^e~ 
3 

4 _ 



_|_ e : 

T* + (e^) 2 + ^ e T*j ^ 



-|- g 1 

By substituting Eq. (40) into Eq. (39), the Eq. (39) can be reduced to: 



(40) 



v(x f t) = 






,V6* -|_g(s* /) 



(41) 



This is the exact solution of the problem, Eq. (33). Table 4 shows the trend of rapid convergence of the 
results of 5 0 (^t) = v^{x } ,t) to S s {x f t) = using the HPM solution toward the exact 
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solution. The maximum relative error of less than 0.037% is achieved in comparison to the exact solution as 
shown in table 4. 



Table 4: 



percentage of relative error (%RE) 






x = .2 


II 


x = .8 


t = 0.1 


s 0 O,t) 


0.1591055626 


0.1376714465 


0.09890297594 




Si(x, t) 


0.007038530707 


0.009754764500 


0.01176933513 




S 2 (x,t) 


0.002482893305 


0.001788283111 


0.0004607861238 




S 3 (x,t) 


0.0001958736131 


0.0002690488417 


0.0002512962487 




s 4 (.x, t) 


0.00004554598541 


0.00002344226507 


0.00001229642845 




S 5 {x,t) 


0.000004998678372 


0.000006399733729 


0.000003681928211 


t = . is 


S a (.x,t) 


0.2155132753 


0.1864801193 


0.1339670591 




S ± (_x, 0 


0.01698529379 


0.02214248709 


0.02558157620 




S 2 (x ! t) 


0.007422698254 


0.005232453593 


0.001127443849 




S 3 (_x,t) 


0.001011675135 


0.001318032849 


0.001182316216 




S 4 (.x, t ) 


0.0003016921205 


0.0001450263132 


0.0001001951821 




■S' 5 ft, 0 


0.00005638430243 


0.00006875985734 


0.00003732434558 


t = .2 


■S' ,3 Or, 0 


0.2605557486 


0.2254546357 


0.1619662983 




■S'lOf, 0 


0.03164334146 


0.03938244479 


0.04388750886 




s 2 (x,t) 


0.01561926168 


0.01076040185 


0.001819030418 




S 3 {x,t) 


0.003225422844 


0.004022798070 


0.003478947590 




s 4 Of, 0 


0.001110362592 


0.0004931283715 


0.0004433744920 




s'bOvO 


0.0003119309814 


0.0003646043520 


0.0001870028440 


Table 


4 shows: 


the percentage 


of relative errors of the 


results of t) = 



5 5 (x t) = Gr, t) of the HPM solution of Eq. (33). 



4. Conclusion: 

In the present work the exact solution of the Newell- Whitehead- Segel nonlinear diffusion equation is 
obtained using the HPM. The validity and effectiveness of the HPM is shown by solving four non-homogenous 
non-linear differential equations and the very rapid convergence to the exact solutions is demonstrated 
numerically. The trend of rapid and monotonic convergence of the solution toward the exact solution is clearly 
shown by obtaining the relative error in compared to the exact solution. The rapid convergence towards the 
exact solutions of HPM indicates that, using the HPM to solve the non-linear differential equations, a reasonable 
less amount of computational work with acceptable accuracy may be sufficient. Moreover it can be concluded 
that the HPM is a very powerful and efficient technique which can construct the exact solution of nonlinear 
differential equations. 



REFERENCES 

Cao, Li. and Bo. Han, 2011. Convergence analysis of the homotopy perturbation method for solving 
nonlinear ill-posed operator equations. Computers and Mathematics with Applications, 61: 2058-2061. 

Ezzati, R. and K. Shakibi, 2011. Using adomian’s decomposition and multi-quadric quasi-interpolation 
methods for solving Newell- Whitehead equation. Procedia Computer Science., 3: 1043-1048. 

Gepreel, A. Khaled, 2011. The homotopy perturbation method applied to the nonlinear fractional 
Kolmogorov-Petrovskii-Piskonov equations. Applied Mathematics Letters, 24: 1428-1434. 

He, J.H., 1998. Approximate analytical solution for seepage flow with fractional derivatives in porous 
media, Comput. Math. Appl. Mech. Eng., 167: 57. 

He, J.H., 1998. Approximate solution of nonlinear differential equations with convolution product 
nonlinearities, Comput. Math. Appl. Mech. Eng., 167: 69. 

1410 



Aust. J. Basic &Appl. Sci., 5(8): 1400-1411, 2011 



He, J.H., 1999. Homotopy perturbation technique. Computer Methods in Applied Mechanics and 
Engineering, 178: 257-262. 

He, J.H., 1999. Int. J. Non-Linear Mech., 344: 699. 

He, J.H., 1999. J. Comput. Math. Appl. Mech. Eng., 17 (8): 257. 

He, J.H., 2000. Int. J. Non-Linear Mech., 351: 37. 

He, J.H., 2005. Application of homotopy perturbation method to nonlinear wave equations. Chaos, 
Solitons & Fractals, 26: 695-700. 

He, J.H., 2005. Application of homotopy perturbation method to nonlinear wave equations. Int. J. 
Nonlinear Sci. Numer.Simul., 6(2): 207. 

He, J.H., 2005. Limit cycle and bifurcation of nonlinear problems, Chaos, Solitons & Fractals. 26 (3): 827. 

He, J.H., 2005. Nonlinear variants of the BBM equation with compact and non-compact physical 
structures, Chaos , Solitons Fractals, 26(3): 695. 

He, J.H., 2005. Phys. Lett. A., 347(4-6): 228. 

He, J.H., 2006. Construction of solitary solution and compaction-like solution by variational iteration 
method, Chaos Solitons Fractals, 29(1): 108. 

He, J.H., 2006. Homotopy perturbation method for solving boundary value problems. Phys. Lett., A 
350:(l-2): 87. 

He, J.H., 2006. Non-perturbative Methods for Strongly Nonlinear Problems, dissertation. de-Verlagim 
IntemetGmbH, Berlin. 

He, J.H., 2006. Some Asymptotic Methods for Strongly Nonlinear Equations Int. J. Mod. Phys., B 20 (10): 
1141. 

Kogak, H., A. Yildirimm, D.H. Zhang and S.T. Mohyud-Din, 2011. The Comparative Boudaker 
Polynomials Expansion Scheme (BPES) and Homotopy Perturbation Method (HPM) for solving a standard non- 
linear second-order boundary value problem. Mathematical and Computer Modelling, 54: 417-422. 

Wazwaz, A.M., 2009. Partial Differential Equations and Solitary Waves Theory. Higher Education Press, 
Beijing and Springer- Verlag Berlin Heidelberg. 

Yildirim, A., 2009. Homotopy perturbation method for the mixed Volterra-Fredholm integral equations. 
Chaos, Solitons& Fractals, 42(5): 2760-2764. 



1411 




